w 


-A18S  820 


THE  RELATIONSHIP  Ot  TENSILE  SPECIMEN  SIZE  AND  GEOMETRV 
EFFECTS  TO  UNIfiUE  <U>  NAVAL  RESEARCH  LAB  WASHINGTON  DC 
P  MATIC  ET  AL  85  OCT  87  NRL-NR-59T6 


UNCLASSIFIED 


F/G  11/6  1 


AD-A186  020 


< 


Naval  Research  Laboratory 


Washington,  DC  20375-5000 


^  H LE  Copy 


NRL  Memorandum  Report  5936 


The  Relationship  of  Tensile  Specimen 
Size  and  Geometry  Effects  to 
Unique  Constitutive  Parameters 
for  Ductile  Materials 


Peter  Matic,  George  C.  Kirby  III  and  Mitchell  I.  Jolles 

Mechanics  of  Materials  Branch 
Material  Science  and  Technology  Division 


October  5,  1987 


DTiC : 


Approved  for  public  release,  distribution  unlimited 


/O' -/•/-/->> /,y«v0 

%  A  -N  J  V  "ft  m%  ,V  a*  _N  -N  -  *  y  n  .  -  _  v  .  •-»*_-“»  -  -  •.■»'«*»*.*#  *v  >  V*  »  ,  /,  < 

fa-f-  a  ill  Val" a  AtJLs  jfl  A!js ifVjfaJLt  iCmJC tkjC ^jCaijC kA  jlK  a/GlA  aAjLI/U  i/^k  LWa  \JV  A  U  Aj 


SECURITY  Classification  of  This  page 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No  0704-0188 


la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2a  SECURITY  CLASSIFICATION  AUTHORITY 


2b  DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 


4  PERFORMING  ORGANIZATION  REPORT  NUM8ER(S) 

NRL  Memorandum  Report  5936 


6a  NAME  OF  PERFORMING  ORGANIZATION  I  6b  OFFICE  SYMBOL  I  7a  NAME  OF  MONITORING  ORGANIZATION 

1  (If  applicable)  j 


lb  RESTRICTIVE  MARKINGS 


3  DISTRIBUTION/AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
unlimited . 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


Naval  Research  Laboratory 


6c  ADDRESS  (City,  State,  and  ZIP  Code) 


Code  6382 


7b  ADDRESS  (City,  State,  and  ZIP  Code) 


Washington,  DC  20375-5000 


8a  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

Office  of  Naval  Research 


8c,  ADDRESS  (City,  State,  and  ZIP  Code) 

Arlington,  VA  22217 


8b  OFFICE  SYMBOL  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(if  applicable ) 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 


61153N22 


PROJECT 

TASK 

WORK  UNIT 

NO  RR022- 

NO 

ACCESSION  NO 

01-48 

DN480-509 

1 1  TITLE  (Include  Security  Classification) 

The  Relationship  of  Tensile  Specimen  Size  and  Geometry  Effects  to  Unique  Constitutive 
Parameters  for  Ductile  Materials 


12  PERSONAL  AUTHOR(S) 

Matic,  Peter,  Kirby,  G.  C.  Ill,  and  Jolles,  Mitchell  I. 


13a  TYPE  OF  REPORT 


13b  TIME  COVERED 


14  DATE  OF  REPORT  (Year,  Month,  Day)  15  PAGE  COUNT 

1987  October  5  41 


COSATI  CODES 


GROUP  SUB-GROUP 


78  SU8JECT  TERMS  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

fconstitutitive 

Tensile  necking 

True  stress-true  strain. 

Fracture  toughness 

Tensile  specimen 

Enerev  density 

19  ABSTRACT  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Straightforward  reduction  of  tensile  specimen  load-displacement  data  to  uniaxial  true  stress-true  strain  data  valid  at 
large  deformations  is  not  possible  for  ductile  material  specimens  which  exhibit  necking.  For  applications  such  as  frac¬ 
ture  mechanics  and  metal  forming  processes,  valid  uniaxial  true  stress-true  strain  data  is  essential  for  the  accurate 
analysis  and  prediction  of  deformations  beyond  those  encountered  at  incipient  necking. 

This  investigation  reports  on  the  use  of  tensile  specimen  size  and  geometry  effects  to  accurately  and  uniquely 
assess  the  uniaxial  true  stress-true  strain  curve.  This  is  done  for  a  ductile  steel  in  the  context  of  an  incremental  elastic 
plastic  constitutive  theory  employing  isotropic  hardening  and  a  von  Mises  yield  function.  Experimental  load-displacement 
data,  experimental  full  field  deformation  data  and  the  corresponding  finite  element  computational  geometries  and  applied 
loads  for  tensile  specimen  length/diameter  ratios  are  known  quantities.  The  uniaxial  true  stress-true  strain  curve  parame¬ 
ters  and  treated  as  unknown  quantities.  Successive  iteration  on  the  solution  curve  resulted  in  successful  correlation 
between  the  experimental  data  and  computational  predictions  for  all  four  specimens  considered.  1  Jt  . 


20  DISTRIBUTION /AVAILABILITY  OF  A8STRACT 
□  UNCLASSIFIED/UNLIMITED  [3  SAME  AS  RPT 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 


Continues 


21  ABSTRACT  SECURITY  CLASSIFICATION 

□  otic  users  UNCLASSIFIED 


22b  TELEPHONE  (Include  Area  Code)  22c  OFFICE  SYMBOL 

(202)  767-5215  Code  6382 


SECURITY  CLASSIFICATION  OF  THIS  PACE 


19.  ABSTRACT  (Continued) 

These  results  support  the  view  that  the  onset  and  evolution  of  tensile  specimen  necking  may  be  accurately 
predicted  by  computational  modeling  of  the  specimen  geometry,  consequence  of  the  interaction  between  these  three 
feature  when  each  is  correcdy  assessed.  This  interaction  was  responsible  for  the  computational  prediction  of  deforma¬ 
tion  symmetry  breaking  in  the  form  of  asymmetric  necking  for  high  length/diameter  ratio  specimens.  This  asymmetry 
was  confirmed  by  the  experimental  data.  As  a  result,  a  priori  assumptions  of  deformation  symmetry  on  the  basis  of 
geometric  symmetry  alone  appear  to  be  incorrect  for  some  geometries. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


V, 


CONTENTS 

INTRODUCTION  . 

EXPERIMENTAL  INVESTIGATION  . 

COMPUTATIONAL  SIMULATIONS  . 

CONTINUUM  MATERIAL  TOUGHNESS  CONCEPTS  . 

SUMMARY  . 

REFERENCES  . 

APPENDIX  I  —  ABAQUS  Finite  Element  Program  Formulation 
APPENDIX  II  —  Laboratory  Specimen  Imperfections  . . 


THE  RELATIONSHIP  OF  TENSILE  SPECIMEN 
SIZE  AND  GEOMETRY  EFFECTS  TO 
UNIQUE  CONSTITUTIVE  PARAMETERS 
FOR  DUCTILE  MATERIALS 

INTRODUCTION 

Accurate  knowledge  of  large  deformation  ductile  material  behavior  is  required  for  many  design 
and  analysis  applications.  Fracture  mechanics  and  metal  forming  processes,  for  example,  fall  into  this 
category.  Two  equally  important  issues  make  the  continuum  characterization  of  inelastic  materials  in 
general  nontrivial.  These  issues  are  (i)  the  rational  development  of  stress,  strain  and  constitutive  for¬ 
mulations  for  a  wide  class  of  materials  and  (ii)  the  determination  of  material  constitutive  parameters 
from  material  test  specimen  response  for  use  in  a  particular  constitutive  formulation.  This  investiga¬ 
tion  will  deal  primarily  with  the  second  issue  as  it  applies  to  a  conventional  incremental  elastic-plastic 
material  formulation  for  large  deformation  and  rotations.  Accomplishing  this  requires  quantitative 
determination  of  a  material’s  uniaxial  continuum  stress  versus  strain  response,  valid  over  its  full  range 
of  deformation  from  zero  to  final  fracture. 

In  the  case  of  ductile  metals,  the  uniaxial  tensile  test  specimen  is  used  to  obtain  global  load- 
displacement  data.  In  principle,  this  data  can  be  easily  normalized  to  uniaxial  material  stress-strain 
data.  From  a  practical  point  of  view,  the  necking  phenomenon  in  ductile  material  specimens  has  been 
a  formidable  barrier  to  accurate  and  complete  constitutive  characterization.  As  a  result,  constitutive 
parameters  are  usually  determined  from  data  obtained  prior  to  necking  over  the  deformation  range 
where  purely  uniaxial  and  homogeneous  deformations  are  assumed  to  exist.  The  resulting  material 
description  is  strictly  valid  only  in  applications  for  which  the  effective  stress  and  strain  measures  do 
not  exceed  the  bounds  of  the  uniaxial  data.  For  ductile  engineering  alloys  presently  in  use.  the  speci¬ 
men  uniaxial  data  obtained  prior  to  necking  only  describes  a  small  portion  of  the  full  strain  range 
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experienced  by  the  material  from  zero  deformation  to  fracture  initiation  at  the  center  of  the  neck. 
This  is  easily  estimated  by  comparing  the  uniaxial  true  strain  at  the  onset  of  necking  with  the  axial 
component  of  strain  at  the  neck  (obtained  from  the  initial  to  final  cross  section  area  ratio)  when  frac¬ 
ture  occurs. 

Traditional  lines  of  investigation  have  attributed  the  onset  and  subsequent  development  of  neck¬ 
ing  to  either  (i)  the  global  instability  of  deformation  predicted  by  one  dimensional  formulations,  (ii) 
the  presence  of  specimen  geometry  imperfections,  usually  in  the  form  of  a  specimen  taper  at  the 
desired  neck  location,  or  (iii)  material  instability  criteria  which  serve  to  localize  deformation.  These 
approaches  to  explain  necking  phenomenology  are  summarized  in  a  review  article  by  Hutchinson 
(1979).  Recent  applications  of  these  approaches  reported  in  the  literature  include  Hutchinson  et.  al. 
(1983)  on  neck  propagation  in  polymers  and  metals,  Tvergaard  et.  al.  (1984)  on  cup-cone  tensile  frac¬ 
ture,  Kleiber  (1986)  on  plastic  localization  in  plain  strain  and  axisymmetric  tension  and  Needleman  et 
al.  (1986)  on  necking  and  failure  in  porous  plastic  solids.  Material  constitutive  parameters  are.  how¬ 
ever,  typically  obtained  from  specimen  data  prior  to  necking  for  use  in  the  analysis  of  necking 
phenomenon.  In  a  practical  sense,  this  defeats  the  purpose  of  the  tensile  specimen  as  a  means  by 
which  to  obtain  constitutive  parameters  valid  over  the  full  range  of  material  deformation. 

Norris  et  al  (1978)  developed  a  uniaxial  true  stress-true  strain  relationship  through  computational 
simulation  of  a  tensile  specimen  gage  length.  Successive  iterations  on  the  true  stress-true  strain  rela¬ 
tionship  were  successful  in  correlating  predicted  specimen  response  with  laboratory  specimen  data. 
The  computational  model  invoked  a  linear  specimen  taper,  however,  to  initiate  necking  at  the  center 
of  the  specimen. 

A  recent  investigation  (Matic,  1985)  suggests  that  the  onset  and  evolution  of  ductile  tensile 
specimen  necking  is  a  natural  consequence  of  the  interaction  between  specimen  geometry,  the  physical 
method  of  load  application  and  material  nonlinearity  No  local  instability,  global  instability  or 
geometric  imperfections  were  used  to  initiate  necking.  The  uniaxial  true  stress-true  strain  constitutive 


parameters  for  the  behavior  of  the  material  at  the  continuum  scale  are  treated  as  the  unknowns  in  the 
problem.  The  solution  curve  was  constructed  to  satisfy  the  requirement  that  computational  simula¬ 
tions  of  the  tensile  specimen  compare  favorably  with  the  measured  behavior  of  the  specimen  in  the 
laboratory.  Examination  of  the  simulation  results  indicated  second  order  nonuniformity  in  the  stress, 
strain  and  energy  fields  which  develop  within  the  specimen  gage  length  immediately  upon  application 
of  the  load.  These  nonuniformities  persist  and  intensify  during  subsequent  elastic  and  inelastic  defor¬ 
mation  of  the  material.  Neck  development  results  after  sufficient  specimen  elongation.  Lateral  con¬ 
traction  and  localization  of  deformation  continue  to  intensify  in  the  vicinity  of  the  neck.  Outside  of 
the  neck  region,  however,  the  material  unloads  as  the  global  reaction  force  decreases  in  response  to 
the  diminishing  cross  sectional  area  at  the  neck. 

The  presence  of  deformation  inhomogeneity  at  the  onset  of  specimen  yielding  is  consistent  with 
earlier  results.  Analytical  results  of  Filon  (1902)  showed  that  boundary  conditions,  representative  of 
physical  load  application,  produce  appreciable  stress  and  strain  inhomogeneties  in  elastic  circular 
cylinders  subjected  to  tension  or  compression  loading.  Needleman  (1972)  reported  inhomogeneous 
plastic  deformation  in  circular  cylinders  at  the  onset  of  yielding.  The  boundary  conditions  in  that 
study  included  radial  constraint  at  the  ends  of  the  cylinder. 

A  subsequent  computational  investigation  (Matic  and  Jolles,  1986)  was  performed  to  examine 
the  sensitivity  of  tensile  specimen  response  to  the  size  and  location  of  a  small  spherical  void  defect  on 
the  axis  of  the  specimen.  The  defect  was  an  order  of  magnitude  smaller  than  the  smallest  relevant 
length  scale  of  the  specimen,  which  in  this  case  was  the  specimen  diameter.  The  results  of  this  inves¬ 
tigation  suggest  that  defects  of  this  type  and  size  can  be  readily  responsible  for  asymmetric  deforma¬ 
tion  and  reduced  apparent  ductility  of  the  specimen. 

The  results  of  these  two  studies  suggested  that  a  hybrid  computational-experimental  technique  is 
a  useful  methodology,  of  sufficient  sensitivity  to  the  physical  parameters  of  the  problem,  with  which 
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to  determine  ductile  material  constitutive  parameters.  The  combined  interaction  between  specimen 
geometry,  applied  load  and  material  behavior,  which  in  combination  produce  specimen  necking,  can 
be  used  to  effectively  solve  for  the  material  constitutive  behavior. 

This  investigation  expanded  upon  these  earlier  results  to  exploit  specimen  size  and  geometry 
effects  in  conjunction  with  the  computational-experimental  approach  to  determine  ductile  material  con¬ 
stitutive  parameters.  In  this  way,  the  relative  sensitivity  of  the  different  tensile  specimen  global 
responses  to  the  local  material  response  could  be  evaluated  and  the  ability  of  a  single,  unique  uniaxial 
true  stress-true  strain  curve  to  predict  the  deformation  of  different  specimens  demonstrated.  Unique¬ 
ness  will  be  referred  to  in  the  sense  that  one  uniaxial  continuum  true  stress-true  strain  curve  can 
predict  the  behavior  of  laboratory  specimens  with  different  geometries.  Twelve  different  tensile 
specimen  geometries  of  HY-100  steel  were  tested  in  the  laboratory.  The  twelve  geometries  were  gen¬ 
erated  from  four  different  gage  lengths  and  three  different  diameters.  From  these  specimens,  the  four 
with  extreme  lengths  and  diameters  were  computationally  simulated  to  develop  a  unique  solution 
curve  for  the  uniaxial  true  stress-true  strain  behavior  which  would  adequately  predict  the  correspond¬ 
ing  observed  specimen  behaviors. 

In  summary,  the  following  approach  was  used: 

(1)  Test  a  family  of  uniaxial  tensile  specimens,  featuring  four  different  lengths  and  three  dif¬ 
ferent  diameters,  to  failure. 

(2)  Treat  the  uniaxial  continuum  true  stress-true  strain  curve  as  the  unknown  quantity  for  com¬ 
putational  simulations  of  the  extreme  specimen  geometries. 

(3)  Generate  a  trial  true  stress-true  strain  curve  for  the  computational  simulations. 

(4)  Compare  the  predicted  axiai  load-displacement  response  and  full  field  lateral  contraction  of 
the  tensile  specimens  with  corresponding  data  from  experiment. 


(5)  Generate  revisions  for  the  true  stress-true  strain  curve  if  differences  between  the  computa¬ 
tional  simulations  and  experiments  are  greater  than  desired  by  repeating  steps  (3)  and  (4). 

(6)  Adopt  the  final  iterate  of  the  true  stress-true  strain  solution  curve  as  the  continuum  behavior 
of  the  material. 

EXPERIMENTAL  INVESTIGATION 

A  family  of  tensile  specimens  was  designed  to  obtain  parametric  data  on  the  effects  of  specimen 
size  and  geometry  for  HY-100  steel.  A  total  of  twelve  different  specimen  geometries  were  examined 
(Fig.  1).  Three  specimens  of  each  geometry  were  tested,  plus  two  extra  specimens  for  the  extreme 
geometries.  Three  different  gage  section  diameters  of  0.76,  1.27  and  1.78  cm  (0.30,  0.50  and 
0.70in)  and  four  different  gage  lengths  of  1.27,  2.54,  3.81  and  5.08  cm  (0.50,  1.00,  1.50  and 
2.00in)  were  used  to  generate  the  twelve  specimen  geometries.  These  geometries  may  be  described 
by  the  gage  length  to  diameter  ratio  (L0/D 0).  A  transition  section  was  included  in  all  specimens 
between  the  gage  section  and  the  2.54  cm  (l.OOin)  diameter  grip  sections  over  which  the  load  was 
applied  at  each  end  of  the  specimen.  The  transition  section  diameter  was  sized  such  that  its  cross  sec¬ 
tional  area  was  double  the  gage  section  area.  For  the  material  under  consideration,  this  ensured  elastic 
response  outside  the  gage  section,  geometrically  similar  transitions  for  all  specimens  and  more  than 
1.0'  axial  distance  between  grip  contact  and  gage  section. 

A  closed-loop  servohydraulic  test  machine  was  used  to  perform  the  tensile  tests.  The  specimen 
ends  were  gripped  by  hydraulic  grips  with  a  grip  length  of  12.3  cm  (3.5in).  All  tests  were  performed 
under  stroke  control  at  a  rate  of  0.127  cm/min  (0.050in/min)  and  a  ramp  loading  function.  Prior  to 
testing,  the  alignment  of  the  grips  was  checked  and  was  found  to  be  within  0.0051  cm  (0.0020in>. 
Thus,  no  significant  bending  was  introduced. 


specimen  diameter,  D0 


0.76  cm  1.27  cm  1.78  cm 


Fig.  1  —  Parametric  Tensile  Specimen  Geometries 


A  clip-on  extensometer  and  a  load  cell  that  is  an  integral  part  of  the  load  frame  were  used  to 
measure  the  elongation  of  the  specimen  in  the  gage  length  and  the  applied  load.  The  calibration  of 
the  load  cell  was  checked  periodically  by  shunt  calibration.  The  extensometer  was  calibrated  directly 
by  a  micrometer  and  a  digital  voltmeter.  The  extensometer  calibration  was  performed  whenever  a 
different  gage  length  was  used.  The  voltage  readings  from  the  extensometer  and  the  load  cell  were 
monitored  by  two  digital  voltmeters  and  converted  to  digital  format.  The  voltmeters  were  simultane¬ 
ously  triggered  to  take  data  by  a  microcomputer  at  1  second  intervals.  The  microcomputer  read  the 
digital  voltages  from  the  voltmeters,  applied  the  appropriate  calibration  factors  and  wrote  the  reading 
number,  load  and  extensometer  displacement  to  a  disk.  The  number  of  data  points  per  test  ranged 


between  200  and  300.  The  microcomputer  was  also  used  to  start  loading  of  the  specimen  by  trigger¬ 
ing  a  function  generator.  On  screen  graphics  allowed  the  test  operator  to  monitor  the  progress  of  the 
test  in  real  time. 

A  photograph  of  the  specimen's  gage  length  was  taken  close  to  the  occurrence  of  the  ultimate 
load.  Subsequent  photographs  were  taken  in  5  to  10  second  intervals  until  failure.  As  failure  became 
imminent,  the  time  interval  between  photographs  was  decreased.  These  photographs  provided  full 
field  specimen  deformation  information  and  a  means  to  characterize  the  neck  region  after  the  test  was 
over.  Since  testing  was  not  halted  as  the  photographs  were  taken,  no  load  relaxation  occurred.  How¬ 
ever,  the  corresponding  load — displacement  data  was  flagged  so  that  the  photograph  could  be  matched 
with  the  proper  specimen  elongation  at  a  later  time.  Typically,  12  photographs  were  taken  for  each 
specimen. 

For  each  set  of  specimens  of  identical  geometry,  the  load  (P)  versus  gage  length  elongation  (5) 
responses  were  essentially  identical  from  the  three  specimens  tested  for  each  of  the  twelve  distinct 
geometries,  demonstrating  data  reproducibility.  The  P  —8  data  was  normalized  to  engineering  stress- 
engineering  strain,  i.e. 


a  = 


P_ 


and 


(1) 


€  = 


(2) 


_5_ 


(31 


where  /t0  is  the  underformed  specimen  cross  section  area  and  L0  is  the  undeformed  specimen  gage 
length.  The  data  is  plotted  in  Fig.  2  for  each  gage  section  diameter  with  gage  section  length  as  a 


parameter.  The  influence  of  specimen  geometry  is  apparent.  The  geometric  similitude  of  the  neck 
geometry  over  all  specimens  (Fig.  3),  with  the  exception  of  the  specimen  with  the  smallest  L0/D ,, 
ratio  of  0.71  is  consistent  with  the  increase  in  global  specimen  ductility  as  the  specimen  gage  section 
is  decreased  for  a  constant  specimen  diameter. 


Fig.  2  -  Engineering  stress-engineering  strain  curves  from  experiment.  Curves  A.  B.  C  and  D  refer  to 
specimen  gage  lengths,  Z.0,  of  1  27,  2.54,  3.81  and  5.08  cm.  respectively  Specimen  diameter.  D,„  of 
(a)  0.76  cm,  (b)  1.27  cm  and  (c)  1.78  cm. 


(a'  (b)  (c)  (d) 

Fig  3  —  Photograph  of  specimen  deformation  at  the  onset  of  speci¬ 
men  failure  for  the  extreme  gage  section  diameters.  Dn.  and 
lengths.  L,,.  producing  aspect  ratios.  l.J Du.  of  (a)  0.71.  (h>  1.67. 

<c>  2.86.  and  < d )  6.67 

(  ( )Ml»l  T  AT  IONA  L  SIMULATIONS 

I'lie  tour  extreme  specimen  geometries,  i.e.  the  specimens  featuring  0.76  cm  or  1.78  cm  (0.3()m 
or  0.70in>  diameter  and  1.27  cm  or  5.08  cm  (0.50in  or  2.00in)  gage  length,  were  selected  for  compu¬ 
tational  simulation.  If  satisfactory  predictions  of  these  four  specimen  responses  could  be  obtained 
using  one  unique  uniaxial  true  stress-true  strain  curve,  then  the  size  and  geometry  independence  ot  the 
material  constitutive  characterization  would  be  demonstrated. 

file  ABAQl'S  finite  element  code  (Hibbitt  et  al.  1984b)  was  used  for  the  computational  simula¬ 
tion  The  finite  element  grids  for  the  four  specimens  are  shown  in  Fig.  4.  Axisymmetry  ot  the 
geometry  required  t hot  only  one-quarter  of  the  specimen  be  modeled. 

Type  C  AX8H  axisymnictric  elements  were  used.  These  are  8-noded  elements  with  quadtalic 
displacement  interpolation  and  an  independently  interpolated  linear  hydrostatic  stress.  The  hydrostatic 
-tress  I-  coupled  to  the  constitutive  relation  using  a  Lagrange  multiplier.  The  use  of  these  hybrid  ele 
merits  prevented  physically  unrealistic  displacement  constraints  from  propagating  through  the  end 
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Such  constraints  can  lead  to  artificial!}  “stiff"  responses  of  standard  elements  for  incompressible  or 
nearly  incompressible  deformations. 


’  'h  •  D,  s)  76  cm  [j,  ’  .’8  ••  D,  -  C76«.rr 

b  . 

l- ii!  4  -  f  inite  element  meshes  for  computation aJ  simulations  of 
specimens  with  aspect  ratios.  Lu/l)u.  of  ta>  0.71.  t h)  1  67.  (cl 
2.X6  and  (d)  6.67 

All  analyses  were  performed  with  full  geometric  nonlinearity  to  account  for  large  strains  and 
large  rotations.  An  updated  l.agrangiati  formulation  is  used  for  incremented  solutions  in  ABACI'S 
Kik's  algorithm,  which  is  load-displacement  controlled,  was  used  to  ensure  numerical  stability,  since 
the  load  decreases  lor  large  specimen  elongations. 

Displacement  boundary  conditions  were  prescribed  to  simulate  the  physical  loading  conditions  of 
the  laboratory  tests.  A  6  J5  cm  (2.50in)  wide  strip  on  the  lateral  surface  of  the  grip  section  ot  each 
specimen  was  subiected  to  uniform  axial  displacements.  These  conditions  reproduced,  lor  computa 
turn,  the  physical  nature  of  the  loads  applied  by  the  hydraulic  grips 

An  incremental  rate  independent  plasticitx  theory  was  used  tor  the  material  constitutive  model 
iHibbilt  et  al  l9K4a).  Total  strains  in  the  multiaxial  strain  state  are  obtained  by  integration  ol  the 


linearly  decomposed  elastic  and  plastic  components  of  the  rate  of  deformation  tensor  DtJ.  (Refer  to 
Appendix  I  for  a  discussion  of  the  ABAQUS  finite  element  program  formulation.)  This  integration  is 
performed  under  the  assumption  that  the  elastic  strains  remain  infinitesimal,  as  is  the  case  in  this 
investigation.  The  total  multiaxial  strain  state  e:j ,  expressed  in  terms  of  elastic  and  plastic  com¬ 
ponents,  is 


+  ifj  ■  (4) 

The  total  logarithmic  uniaxial  strain  t ,  consistent  with  the  integration  of  the  rate  of  deformation  tensor 

for  a  multiaxial  strain  state,  is  decomposed  as 


e  =  ef  +  6P 


The  yield  function  /  takes  the  form 


f(Tij)  -  TfeP)  (6) 

where  r(;  and  t  are  multiaxial  and  uniaxial  Kirchoff  (or  Treffetz)  stress  states,  respectively.  The 
associated  flow  rule  governs  plastic  strain  increments  by  the  relation 


dt' = x£- 


In  the  case  of  purely  elastic  behavior  X  =  0.  For  active  yielding. 


X  >  0 


X  =  dep  -7- - - 


Plastic  strain  increments  also  satisfy  a  dissipation  equivalence  condition 


t  dep  =  t,j  dep 


1 1 


and  a  consistency  condition 


kv.' 


4*-  dr,.  -  —  dfP  =  0.  (11) 

Bt.j  3fP 

The  von  Mises  yield  function 

f(T,j)  =  ~  ( SjjSjj )l/2  (12) 

was  employed.  The  deviatoric  stress  tensor  s,j  is  defined  as 

Sjj  —  Tjj  ~  —  T/Jc&ij  .  (13) 

where  the  hydrostatic  component  of  stress  is  rtt/3. 

The  Kirchoff  stress  and  logarithmic  strain  measures  are  employed  because  of  the  advantages 
gained  in  computational  implementation.  The  Kirchoff  stress  tensor  t,j  is  approximately  equal  to  the 
more  physically  motivated  Cauchy  stress  tensor  <ri;  for  deformations  involving  only  small  changes  in 
volume.  This  condition  is  satisfied  in  these  analyses.  The  uniaxial  true  stress-true  strain  constitutive 
response  of  the  material  is  formally  input  as  Cauchy  stress  and  logarithmic  strain  for  the  ABAQUS 
program. 

The  form  of  the  uniaxial  true  stress-true  strain  curve  was  subject  only  to  the  restrictions  that  it 
be  monotonically  increasing  and  convex.  These  restrictions  are  consistent  with  Drucker's  postulates 
for  elastic-plastic  constitutive  formulations  of  the  type  employed  here.  The  curve  itself  is  expressed 
in  a  multilinear  form  as  stress  and  plastic  strain  pairs  supplementing  the  elastic  modulus  and  Poisson's 
ratio  for  the  material. 

The  postulated  material  true  stress-true  strain  curve  was  evaluated  by  one  or  two  of  the  finite 
element  models  per  iteration  in  order  to  assess  its  ability  to  satisfactorily  predict  the  observed 
responses  of  these  tensile  specimens.  The  following  guidelines  were  developed  from  the  computa¬ 
tional  simulations  conducted  in  this  investigation: 
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(1)  Low  L0/D o  specimens  provide  a  good  gross  indication  of  the  material  true  stress-true  strain 
curve.  A  larger  relative  portion  of  the  specimen  is  part  of  the  necking  process,  and  as  a  result  is 
under  monotonically  increasing  deformation. 

(2)  High  L0/D0  specimens  are  most  useful  to  guide  minor  adjustment  to  the  curve  obtained  from 
the  low  Lq/Dq  specimen  data.  A  smaller  portion  of  the  higher  Lq/Dq  specimens  is  involved  in  the 
actively  yielding  neck  in  the  specimen.  As  a  result,  the  balance  between  the  values  of  stress,  strain 
and  tangent  modulus  of  the  curve  and  the  relative  tendency  for  the  neck  to  continue  deforming  while 
the  remainder  of  the  specimen  unloads  are  addressed. 

(3)  Translation  of  the  material  true  stress-true  strain  curve  in  the  stress  axis  direction  results  in  a 
comparable  translation  of  the  predicted  engineering  stress-engineering  strain  curve. 

(4)  Expansion  of  the  material  true  stress-true  strain  curve  in  the  strain  axis  direction, 
corresponding  to  deformation  prior  to  peak  engineering  stress,  tends  to  accelerate  the  onset  of  neck¬ 
ing. 

(5)  Expansion  of  the  material  true  stress-true  strain  curve  in  the  strain  direction,  corresponding 
to  deformation  subsequent  to  peak  engineering  stress,  tends  to  accelerate  the  advanced  stages  of  neck¬ 
ing. 

After  approximately  five  iterations  the  material  true  stress-true  strain  curve,  developed  using  the 
guidelines  listed  above,  was  applied  to  a  full  series  of  analyses  performed  for  all  four  specimen 
geometries  modeled.  The  predicted  engineering  stress-engineering  strain  curves  were  in  good  agree¬ 
ment  for  the  three  specimens  with  the  lowest  L0/D0  ratios  of  0.71.  1.67  and  2.33  from  the  four 
extreme  geometries  as  used  in  the  development  of  the  material  true  stress-true  strain  curve  (Fig.  5). 
The  specimen  with  the  highest  L0/D0  ratio  of  6.67  predicted  higher  engineering  strains  than  antici¬ 
pated  based  upon  the  half-gage  length  elongations  provided  by  the  analyses.  In  fact,  the  predicted 


engineering  stress-engineering  strain  curve  for  the  high  L0/D0  specimen  was  close  to  that  for  the 
second  highest  L0/D0  specimen.  This  apparent  anomaly  of  one  predicted  specimen  response  out  of 
the  four,  in  light  of  the  observed  relative  sensitivity  of  the  predicted  specimen  responses  to  the 
material  true  stress-true  strain  curve,  was  contrary  to  what  would  be  anticipated. 


Fig.  5  —  Engineering  stress-engineering  strain  curves 
measured  from  experiment  and  predicted  by  computa¬ 
tional  simulation  allowing  for  symmetric  deformation 
only:  0,  LJD0  -  0.71;  +  ,  L</D0  -  1.67;  □,  LJ Da 
-  2.86;  X,  L<JD0  -  6.67. 

Upon  examination  of  deformed  shapes  (Fig.  6)  of  the  specimens  from  the  analysis,  it  became 
apparent  that  symmetric  double  necking  of  the  high  L0/D0  ratio  specimen  was  being  predicted,  while 
all  three  low  L0/D0  specimen  analyses  were  predicting  symmetric  single  necking.  This  translated 
directly  into  a  higher  engineering  strain  for  the  high  L0/D0  specimen  since  one  whole  neck  was 
present  in  the  half-gage  length  model,  as  compared  to  one-half  neck  per  half-gage  length  for  the  three 
specimens  with  lower  L0/D0  ratios.  A  neck  region  is,  of  course,  the  source  of  particularly  high  axial 
elongations. 
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Fig.  6  —  Predicted  deformation  at  the  onset  of  specimen  failure  for 
aspect  ratio,  L^D0,  of  (a)  0.71,  (b)  1.67,  (c)  2.86  and  (d)  6.67 

Symmetric  double  necking  is  not  physically  observed  in  the  laboratory.  For  it  to  actually  occur, 
both  spatial  homogeneity  and  temporal  simultaneity  of  an  extremely  high  order  would  have  to  be 
present.  Without  such  homogeneity  and  simultaneity,  a  single  neck  would  develop  asymmetrically,  as 
one  necking  site  would  be  preferred  over  another,  resulting  in  monotonically  increased  loading  at  the 
actual  necking  site  and  load  reduction  at  the  other  site,  consistent  with  the  results  of  the  other  speci¬ 
mens  in  this  study.  In  order  to  translate  the  symmetric  double  neck  results  into  asymmetric  single 
neck  results,  the  influence  of  one  neck  was  removed  in  the  following  manner. 

Surface  stretch  ratios  were  approximated  using  initial  nodal  coordinates  and  nodal  displacements 
at  the  surface  (Fig.  7).  At  the  neck  locations  the  strains  were  significantly  higher  than  those  along  the 
remainder  of  the  specimens.  To  obtain  an  approximation  to  the  total  asymmetric  specimen  length,  the 
contribution  of  strains  at  one  neck  to  the  total  specimen  length  were  subtracted  from  the  total.  The 
average  “background”  strain  contributions  to  the  specimen  length  were  substituted  for  the  neck  seg¬ 
ment  just  removed.  The  result  was  a  total  specimen  length  based  on  the  contributions  of  one  neck 
region  and  a  complementary  background  deformation.  Thus, 


where  L,„  is  the  length  of  the  specimen  with  a  single  neck,  L ^  is  the  length  of  the  specimen  with  a 
double  neck.  Ln  is  the  length  contribution  of  a  single  neck  and  Lb  is  the  length  contribution  due  to 
background  deformation  outside  the  neck.  The  engineering  strain  for  the  specimen  will  be 


This  resulted  in  good  agreement  between  the  engineering  stress-engineering  strain  curve  predicted  by 
the  analysis  and  the  experimental  data  for  the  high  L0/D0  geometry  as  seen  in  Fig.  8. 


Fig.  7  —  Translation  of  predicted  symmetric  double  neck  specimen 
response  into  asymmetric  single  neck  specimen  response 


For  this  interpretation  to  be  correct,  the  laboratory  tensile  specimens  should  show  evidence  of 
symmetric  deformation  for  the  three  lower  L  ID  specimens  and  asymmetric  deformation  for  the  high 
L/D  specimen.  An  eccentricity  measure  to  quantitatively  describe  necking  asymmetry  may  be 
defined  by  the  expression 
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Fig.  8  —  Engineering  stress-engineering  strain  curves  measured 
from  experiment  and  predicted  by  computational  simulation 
allowing  for  symmetric  or  asymmetric  specimen  deformation: 
O.  L0/D0  =  0.71;  +.  L  0/D0  =  1.67;  □.  L0/Da  =  2.86;  X. 
Lq/Dq  ~  6.67. 


e  = 


(16) 


where  e  is  the  eccentricity,  L  is  the  length  of  the  specimen  gage  length  just  prior  to  specimen  fracture 
and  z„  is  the  distance  from  a  specimen  gage  mark  to  the  point  of  minimum  specimen  diameter. 
Clearly,  when  the  neck  occurs  at  the  center  of  the  specimen,  z„  -  L/ 2  and  therefore  e  -  0.  When 
the  neck  is  at  either  end  of  the  gage  length,  z„  =  0  or  L  and  therefore  e  =  1.  Perfect  symmetry 
coincides  with  e  =  0,  while  maximum  possible  asymmetry  coincides  with  e  =  1 .  These  results  are 
plotted  in  Fig.  9.  For  Lq/Dq  less  than  1.67,  the  experimental  data  shows  very  low  values  of  e .  indi¬ 
cating  almost  perfect  symmetry.  For  Lq/Dq  greater  than  1.67,  the  eccentricity  values  tend  to  increase 
to  an  average  of  about  0,5  for  Lq/Dq  greater  than  4.0. 


The  computational  prediction  of  asymmetry  in  specimen  deformation  is  verified  by  the  experi¬ 
mental  data.  The  computational  predictions,  performed  for  ideal  specimens  free  of  imperfections, 
establish  eccentricity  values  for  comparison  purposes.  The  maximum  measured  diameter  imperfec- 


17 


•TV f 


□  Experimental  Data 
O  Computational  Prediction 


0.61 


0.4(-  g  0 


□ 


□ 

o 


1 


o 

□ 

□ 

o 


Q 


_1_ 

6.0 


_1 

8.0 


speciman  aspect  ratio,  Lq/Dq 


Fig.  9  —  Eccentricity,  e,  of  specimen  deformation  at  failure 
versus  specimen  aspect  ratio,  LJDq.  □,  data  from  experi¬ 
ments;  O,  predictions  from  computational  simulations. 


tions  exhibited  by  the  laboratory  specimens,  tabulated  in  Appendix  II,  were  on  the  order  of  10~2. 
Imperfections  of  this  magnitude,  in  line  with  the  conclusions  reached  by  Matic  and  Jolles  (1986)  on 
the  role  of  geometric  imperfections  in  asymmetric  necking,  are  consistent  with  the  measured  neck 
eccentricity  values  and  scatter  of  the  laboratory  data  as  compared  to  the  computationally  predicted 
values. 

The  hierarchy  of  imperfections  effects  is  apparent: 

(1)  Imperfections  are  not  necessary  for  necking  to  occur 

(2)  Imperfections  which  are  effectively  one  order  of  magnitude  smaller  than  the  effective  speci¬ 
men  dimension  (in  this  case  diameter)  are  responsible  for  asymmetric  necking 

(3)  Symmetry  imperfections  significantly  smaller  than  the  effective  imperfection  size  are  suffi¬ 


cient  to  determine  precedence  between  two  apparently  simultaneous  events  at  the  continuum  scale 


The  solution  curve  is  shown  in  Fig.  10.  The  tabular  form  for  the  curve  is  given  in  Table  1. 
The  uniaxial  continuum  true  stress-true  strain  curve  show  a  uniaxial  yield  stress  of  7.17  x  108  N/m: 
(104,000  lb/in~)  and  a  uniaxial  maximum  stress  of  just  under  1.03  x  109  N/m:  (150,000  lb/in;).  The 
negligible  slope  at  the  higher  strain  value  is  consistent  with  the  ductility  of  HY-100.  The  notion  of  a 
limit  load  in  a  rigid — perfectly  plastic  sense  can  be  viewed  as  extrapolation  back  to  the  stress  axis 
from  the  maximum  stress  value. 
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Fig.  10  —  Solution  curve  for  uniaxial  continuum  true  stress-true 
strain  response  of  HY-100  steel 


The  photographs  taken  of  the  tensile  specimens  strongly  suggest  the  concept  of  a  fully  developed 
neck.  Geometric  similitude  is  in  evidence  for  all  but  the  specimens  with  a  gage  length  of  1.27  cm 
(0.50in)  and  a  diameter  of  1.78  cm  (0.70in).  This  specimen  geometry  had  the  lowest  L0/DQ  ratio, 
and  featured  the  lowest  percent  reduction  of  area  on  the  average.  Further  examination  of  this 
geometric  similitude  in  the  neck  profile  tends  to  support  a  view  that  the  length  of  the  neck  remains  a 
constant  multiple  of  the  original  specimen  diameter.  The  inflection  points  of  the  deformed  specimen 
contour,  which  appear  above  and  below  the  minimum  specimen  diameter  at  the  neck,  remain 
separated  by  a  constant  distance  A  from  incipient  necking  to  specimen  fracture  (Fig.  1 1). 
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The  inflection  points,  while  providing  the  intrinsic  geometric  reference  length  for  the  neck,  do 
not  fully  distinguish  the  actively  flowing  material  of  the  neck  from  the  material  outside  from  the  neck. 
Outside  the  neck,  flow  has  ceased  and  unloading  is  in  progress  as  the  axial  stress  acting  on  the  speci¬ 
men  cross  section  decreases  with  the  global  load  due  to  neck  development.  In  fact,  for  the  HY-100 
I  specimens  of  L0/D0  ratio  greater  than  one,  additional  reference  lengths  based  on  the  stress  (or  strain) 

states  in  the  neck  provide  a  more  complete  picture  on  which  to  define  a  neck  length.  Specifically,  the 
hydrostatic  stress  contours  which  are  bounded  and  contain  the  center  of  the  neck  extend  over  a  total 
length  2A  centered  on  the  minimum  diameter  plane.  The  von  Mises  equivalent  stress,  which  under 
the  assumptions  of  plasticity  theory  as  employed  in  the  computational  simulations  govern  the  plastic 
deformation,  is  essentially  uniform  over  a  distance  3A.  centered  on  the  plane  of  minimum  diameter. 

i 

The  neck  radius  at  both  ends  of  this  distance  is  almost  that  of  the  specimen  outside  the  neck. 

From  a  phenomenological  point  of  view,  this  last  observation  provides  a  functional  definition  of 
|  the  neck  based  on  an  intrinsic  length  over  which  energy  is  actively  dissipated.  For  the  particular 

material  under  consideration  this  length  is  3A,  or  three  times  greater  than  the  intrinsic  geometric 
length  parameter  A  defined  from  the  inflection  point  separation  on  the  specimen  profile.  Thus,  prior 

f 

i  to  the  onset  of  necking  the  energy  dissipation  throughout  the  specimen  is  driven  by  axial  extension 

As  the  neck  develops  energy  dissipation  is  confined  to  the  neck,  promoted  by  lateral  contraction  of 
the  neck. 

Working  backward  from  the  final  specimen  profile,  the  initial  length  of  the  neck  on  the 

undeformed  specimen  can  be  calculated,  and  the  ratio  (L0)neck/£>0  determined.  For  the  HY-100 
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material  under  consideration. 


This  implies  that  tor  a  fully  developed  neck  to  form. 


„8) 

D0  D  o 

That  is  to  say,  the  specimen  aspect  ratio  must  be  at  least  that  required  for  a  fully  developed  neck. 
The  specimen  L0/D0  ratio  fell  below  1.0  only  for  the  L0  =  1.27  cm  (0.50in).  D0  =  1.78  cm 
(0.70in)  specimens  from  the  full  set  of  over  three  dozen  specimens  featuring  twelve  different 


CONTINUUM  MATERIAL  TOUGHNESS  CONCEPTS 


The  strain  energy  per  unit  mass  at  a  given  instant  during  deformation,  will  be 


1  AW 

w  =  lim - 

av-o  p  Ak 


a  a 

w  =  f  -J-  d(,.  (20) 

Jo  p  ,J 

where  p  is  the  mass  density,  W  is  energy  and  V  is  volume.  This  energy  density  incorporates  both 
stress  and  strain  into  a  fundamental  quantity  relevant  to  thermodynamic  descriptions  of  material  defor¬ 
mation  and  damage.  The  energy  density  is  a  scalar  quantity  which  takes  into  account  all  components 
of  the  stress  and  strain  tensors  in  a  physically  consistent  manner. 


Failure  of  the  material  due  to  fracture,  at  the  continuum  scale,  can  be  associated  with  the  value  of  the 
energy  density  at  which  fracture  occurs.  Thus,  the  material  toughness  is  defined  as 


Vf  =  J  '  ^ 

J0  a 


where  is  the  critical  strain  energy  density  value. 


For  ductile  metals,  the  material  densitv  varies  onlv  slmhtlv,  even  over  large  deformations.  For 


this  reason,  it  is  common  to  define  an  energy  per  unit  volume  density 


w  =  lim 
.11''  —  o 


AW 

AV 


(22) 


=  fo'  a,y  deu  .  (23) 

with  an  associated  critical  value 

^  *  i0  °'J  de‘J  •  (24) 

The  energy  per  unit  mass  is  fundamental,  but  the  energy  per  unit  volume  is  equally  appropriate  for 
constant  volume  deformation  processes. 


The  use  of  an  energy  density  criterion  for  continuum  fracture  is  consistent  with  the  ideas 
expressed  by  Freudenthal  (1950)  on  material  behavior  and  scaling  considerations.  Gillemot  (1976) 
examined  cylinderical  tensile  specimens  using  analytical  and  empirical  techniques.  The  strain  energy 
density  per  unit  volume,  absorbed  by  the  material  up  to  the  instant  of  fracture,  was  calculated  from 
global  specimen  response  and  deformation  geometry.  Calculations  of  critical  strain  energy  density 
values  from  computational  simulations,  however,  are  able  to  address  the  large  deformation  continuum 
toughness  in  detail  appropriate  for  structural  integrity  analysis.  Applications  to  structural  integrity 
prediction,  using  continuum  material  toughness  concepts  as  outlined  here,  include  investigations  by 
Matic  and  Jolles  (1987a,  1987b)  into  the  defect  tolerance  of  welded  components  and  the  role  of 
material,  geometry  and  applied  load  on  weld  system  performance. 

For  a  multiaxial  state  of  stress,  each  of  the  six  stress-strain  pairs,  three  normal  and  three  shear, 
must  be  evaluated  and  summed,  i.e. 

r(e,|)'  J  r (<:A  ,  ,<€iA 

Wc  -  J0  ffll  “£ll  +  j0  a22  dt22  +  )0  CT33  ^f33  + 

+  )0  ff|2de,2  +  fo  <7;3^£:3  +  J0  <ri|<*e3l-  (25) 

It  should  be  noted  that  one  or  more  individual  terms  in  the  multiaxial  expression  can  be  negative. 
Their  total,  wr,  must  be  positive,  however. 
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Tensile  specimen  failure  is  known  to  initiate  at  the  center  of  the  specimen  neck,  where  the 
stresses  and  strains  are  multiaxial  in  character.  The  stresses  and  strains  are  plotted  for  the  two  1  78 
cm  (0.70in)  diameter  specimens  for  finite  element  integration  point  nearest  to  the  center  of  the  speci¬ 
men  (Fig.  12).  Symmetry  conditions  at  the  center  of  these  specimens  are  reflected  in  the  absence  of 
any  significant  shear  stress  or  strain  components.  The  I0/D0  ratios  for  these  specimens  are  0.72  and 
2.86,  respectively.  The  lower  ratio  specimen  is  below  the  threshold  ratio  of  1.0  for  fully  developed 
necking.  The  higher  ratio  specimen  is  above  the  threshold. 

All  strain  magnitudes  are  monotonically  increasing  over  the  deformation  history.  For  the  low 
L0/D0  specimen,  the  absence  of  a  fully  developed  neck  results  in  monotonically  increasing  values  of 
stress  over  the  entire  deformation.  Integration  of  these  curves  yields  a  critical  energy  density  at  the 
fracture  initiation  site  of  8.89  x  108  N-m/m3  (1.29  x  105  lb-in/n3).  For  the  high  LQ/D 0  specimen, 
the  presence  of  a  fully  developed  neck  is  preceded  by  essentially  uniaxial  deformation.  As  a  result, 
the  z  component  of  stress  monotonically  increases  while  the  r  and  d  components  of  stress  are  initially 
zero  and  monotonically  increase  only  after  necking  has  begun.  Integration  yields  a  critical  energy 
density  of  1.20  x  109  N-m/m3  (1.74  x  104  lb-in/n3). 

Essentially  proportional  loading  is  present  in  the  low  L0/D0  specimen  stress  state.  Distinctly 
nonproportional  loading  is  evident  for  the  high  L0/D0  specimen.  The  difference  in  the  critical  energy 
density  is  significant,  particularly  in  light  of  the  relative  insensitivity  of  the  stress-strain  components 
to  different  loading  histories.  The  difference  in  the  two  sets  of  stress-strain  plots  occurs  for  the  r  and 
6  components  and  is  probably  an  indication  of  different  energy  dissipation  mechanisms  being 
developed  during  the  initial  deformation  and  yielding  of  each  specimen.  While  on  the  larger  scale  the 
difference  between  the  plots  may  appear  small,  in  an  energy  sense  it  is  much  greater  The  sensitivity 
of  the  continuum  material  toughness  to  nonproportional  deformation  histories  has  significant  implica¬ 
tions  for  fracture  initiation,  crack  growth  and  structural  integrity  predictions 


SUMMARY 


A  hybrid  computational-experimental  approach  for  the  determination  of  material  uniaxial  consti¬ 
tutive  parameters  is  applied  to  the  interpretation  of  tensile  test  specimens  of  different  sizes  and 
geometries.  The  iterative  development  of  a  single  unique  material  true  stress-true  strain  curve  for 
which  computational  predictions  agree  with  experiment  data  for  different  specimens  was  discussed. 
The  unique  solution  curve  for  HY-100  steel  is  presented. 

The  tensile  specimen  neck  is  observed  to  be  geometrically  invariant  to  size  and  geometry  effects 
if  the  specimen  L0/D0  aspect  ratio  is  greater  than  one.  This  ratio  defines  a  lower  bound  for  a  fully 
developed  neck.  An  upper  bound  for  the  L0/D0  ratio  exists  in  terms  of  specimen  deformation  sym¬ 
metry.  Evidence  for  this  upper  bound  was  observed  computationally  and  confirmed  from  experiment 
data. 

Continuum  material  toughness  concepts  were  discussed  in  terms  of  the  critical  energy  density  to 
produce  fracture.  Computational  reconstruction  of  the  stress-strain  histories  leading  to  fracture  initia¬ 
tion  for  the  HY-100  specimens  strongly  suggest  that  the  critical  energy  density  is  different  for  an  arbi¬ 
trary  nonproportional  deformation  history  when  compared  to  a  proportional  deformation  history.  This 
seems  to  be  true  despite  the  general  insensitivity  of  the  stress-strain  components  to  the  different  load¬ 
ing  histories  encountered. 
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APPENDIX  I 


ABAQUS  Finite  Element  Program  Formulation 

A  brief  summary  of  the  ABAQUS  finite  element  program  is  included  here  for  reference  and 
completeness  in  view  of  the  various  possible  theoretical  formulations  available  for  geometrically  and 
materially  nonlinear  analysis.  The  reader  is  referred  to  Hibbitt  et  al  (1984a)  for  further  details. 

For  a  material  particle  at  reference  position  X,  at  time  t  equal  to  zero,  its  current  position  x,  at  a 
later  time  may  be  expressed  as 


x{  -  Xj(Xi,t).  (Al) 

An  orthogonal  coordinate  system  of  unit  base  vectors  is  assumed.  The  relative  position  dx,  between 
two  material  particles  in  their  current  configuration  may  be  expressed  in  terms  of  their  relative  posi¬ 
tion  in  their  reference  configuration 


dx, 

*  -  W"' 


(A2) 


-  F;  dX> 


(A3) 


where  F,}  is  defined  as  the  deformation  gradient. 


The  velocity  of  the  material  particle  will  be 


V'r 


dt 


(A4) 
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The  incremental  velocity  between  two  particles  will  be 


<3v, 

Jv'  - 


3v,  dXj 


dxj  dXk 


=  L,jFjkdXk 


where  LX]  is  the  velocity  gradient  at  a  material  point  in  the  body.  The  velocity  gradient  may  be 
decomposed  into  the  sum  of  its  symmetric  and  antisymmetric  tensor  components  as 


+vy.,)  +  T(Vi.;  -Vj.i) 


=  A; 


(A10) 


where  is  the  rate  of  deformation  tensor,  associated  with  material  straining,  and  is  the  material 
spin  tensor  associated  with  pure  rotation  at  the  material  point. 

For  a  material  body  of  volume  V  and  surface  S  in  its  current  configuration,  force  equilibrium  in 
integral  form  is  given  by 


Jsf,d5  =  0 


(All) 


where  t,  is  the  surface  traction  vector  per  unit  of  current  area  and  /,  is  the  body  force  vector  per  unit 
of  current  volume.  The  Cauchy  stress  tensor  a,,  is  related  to  t ,  by 


r,  <Ji/nj 


(A  12) 


V  A*" 


where  n;  is  the  outward  unit  vector  normal  to  surface  S.  Using  (A12)  in  (All)  along  with  Gauss' 
theorem  to  transform  the  surface  integral  to  a  volume  integral. 


\v(at]J  +f,)dV  =0.  (A  13) 

The  integral  must,  in  general,  be  zero.  Therefore, 

°ij.j  +/,  =  0.  (A  14) 

Moment  equilibrium  in  integral  form  requires 

$s*ilje‘jkdS  +  \vxJje>jkdV  =  0  ,A15> 

where  e,jk  is  the  permutation  symbol.  Using  Gauss'  theorem 

+  xiajUeijk  +  xi  ajleijk  .1  +  xif j  eijk  )dV  =  0  <A16) 

Evaluating  each  term  in  this  expression  and  using  (A14)  provides  that  <riy  must  be  symmetric,  i.e. 

a,j  =  Oj,  (A  17) 

The  differential  force  equilibrium  expression  (A  14)  is  used  to  generate  a  virtual  work  expression 
by  taking  the  dot  product  with  a  virtual  velocity,  v",  and  integrating  over  the  volume.  This  takes  the 
form,  after  application  of  Gauss'  theorem, 

\v  v-’,v  +/><’]  dV  =  °-  (A18) 

The  virtual  velocity  field  must  satisfy  continuity  and  prescribed  velocity  requirements.  Expressing  the 
first  integral  in  terms  of  tractions  and  rearranging 

\voty,.,dv  =  \st,v;ds +  \j,v;jv  ,ai9) 

From  the  definition  of  the  velocity  grrdient.  the  virtual  velocity  gradient  is 


(A20) 


=  £>,’+«,’  (A21) 

where  the  virtual  rate  of  deformation.  D’.  is  symmetric  and  the  virtual  spin,  Q,’,  is  antisymmetric. 
In  light  of  the  Cauchy  stress  tensor  symmetry,  it  follows  that 

\y°iiDiidV  =  Js'«V*S  +  \yf,v'dV .  (All) 

The  left  hand  side  of  this  expression  is  converted  to  an  integral  taken  over  the  reference  volume  V0 , 
i.e. 

Jsuv  V0ai]D'JJdV0  =  \st,v'dS  +\vflv'dV .  (A23) 

where  J  is  the  determinant  of  the  Jacobian.  J  is  equal  to  the  ratio  of  the  differential  volumes  in  the 
reference  and  current  configuration. 

The  Ki rchoff  (or  Treffetz)  stress  is  defined  as 

r,j  =  J  a ij .  (A24) 

Therefore, 

\vr‘JD'idVo  =  v'dS  +  \j>v’dV  (A25) 

The  Kirchoff  stress,  r,y,  is  equal  to  the  Cauchy  stress.  at],  when  J  is  equal  to  unity.  For  the 
material  behavior  considered  in  this  analysis  J  deviates  from  unity  by  only  a  small  amount.  There¬ 
fore,  the  constitutive  formulation  in  terms  of  Cauchy  stress,  as  described  in  the  main  body  of  this 
paper  will  be  applicable  when  expressed  in  terms  of  Kirchoff  stress. 

Equation  (A26)  is  the  virtual  velocity  equilibrium  equation  and  is  the  basis  for  interpolation  and 
discretization  for  the  finite  element  formulation.  The  virtual  velocity  at  a  point  inside  an  element  will 
be  interpolated  as 
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(A26) 


where  the  summation  is  understood  to  involve  only  those  nodal  values  vj  relevant  to  a  particular  ele¬ 
ment.  The  virtual  rate  of  deformation  becomes 


°ij  =  y  (N,k  vk ) v  +  (Njk  vt’; )j  (Ml) 

=  y  (ty*  j  +  ty* )Vfc  ( A28) 

Therefore,  (A25)  becomes 

\v\  TijWk.j  +  NjkJ)vkdV0  =  J^JV^iS  +  lyfiNafidV  (A29) 

or,  since  vt*  appears  in  each  integral, 

\v  \ t^j  +  NjkJ)dV o  =  ^tjNfidS  +  \vflNlkdV .  (A30) 


The  integrations  are  taken  with  respect  to  the  discretized  volumes  and  surfaces  consistent  with  element 
and  nodal  definitions. 

Equation  (A30)  forms  the  basis  for  the  finite  element  solution.  An  incremental  solution  is  gen¬ 
erated  by  the  Newton  algorithm.  The  Jacobian  of  (A30)  is  required  for  the  incremental  solution,  and 
results  in  the  introduction  of  a  Kirchoff  stress  rate,  r,y,  into  the  incremental  equilibrium  equation. 
The  stress  rate  takes  the  form 

Tij  =  h,  +fl.***;  +  T,k  tljk  (A31) 

where  r,;  is  the  total  stress  rate,  h,  's  the  material  (or  Jaumann)  stress  rate  associated  with  the 
material  response  through  the  constitutive  formulation.  The  remaining  terms  involving  fll;  are  stress 
rate  contributions  due  to  rotation.  In  practice,  all  constitutive  calculations  are  performed  using  r 


The  total  strain  measure  t:j  is  calculated  in  a  manner  consistent  with  the  integration  of  DtJ 
This  generates  logarithmic  strain  measures  which  are  energy  conjugate  to  the  Kirchoff  stress  r, 
which,  as  discussed  above,  is  approximately  equal  to  the  Cauchy  stress  atJ  in  .he  analyses  under  con 


sideration. 


APPENDIX  II 


Laboratory  Specimen  Imperfections 

Prior  to  testing,  the  tensile  test  specimen  gage  section  diameters  were  measured.  Deviations 
from  the  nominal  specimen  gage  diameter  D0  are  tabulated  below  for  each  specimen  as  the  maximum 
and  minimum  specimen  diameters.  (DQ)mAX  and  (D0)mm.  measured  along  the  gage  length.  The  diame¬ 
ter  imperfection  d0  is  defined  as 


(Dp)  max 
(D0) 

max  +  (^o)min 


(AIM) 


This  may  be  interpreted  as  the  average  deviation  (from  the  average  diameter)  divided  by  the  average 


diameter. 


The  specimens  modeled  for  the  computational  simulations  were  assumed  to  be  free  of  any 
geometric  imperfections.  The  scatter  of  neck  eccentricity  values  with  respect  to  the  computational 
predictions  of  neck  eccentricity  values,  as  discussed  above,  are  consistent  with  the  small  imperfections 
in  the  specimen  diameters  listed  in  Table  All- 1 . 


Table  A 1 1  -  i  Maximum  measured  diameier  imperfection  value  of  laboratory  specimens. 
Specimens  not  measured  prior  to  testing  are  denoted  by  “  —  "  values. 


Nominal  Gage 

Nominal  Gage 

Maximum  Gage 

Minimum  Gage 

it 

Length.  1 0 

Diameter.  £),, 

Diameter,  i  D  .y  ^a< 

Diameter.  <£>,,)  min  | 

1.27 

cm 

(0.50 

in) 

0.76 

cm 

(0.30 

in) 

07~0 

cm 

(0.303 

in) 

0765 

cm 

(0.301 

in) 

0.767 

cm 

(0.302 

in) 

0.762 

cm 

(0.300 

in) 

0-772 

cm 

(0.304 

in) 

0767 

cm 

(0.302 

in) 

1.27 

cm 

(0.50 

in) 

1.275 

cm 

(0.502 

in) 

1.270 

cm 

(0.500 

in) 

1.285 

cm 

(0.506 

in) 

1.280 

cm 

(0.504 

in) 

1.275 

cm 

(0.502 

in) 

1.275 

cm 

(0.502 

in) 

1.78 

cm 

(0.70 

in) 

1.778 

cm 

<0.700 

in) 

1.778 

cm 

(0700 

in) 

1.791 

cm 

(0.705 

in) 

1786 

cm 

(0703 

in) 

1.786 

cm 

(0.703 

in) 

1.781 

cm 

(0.701 

in) 

1.773 

cm 

(0  698 

in) 

1.770 

cm 

(0.697 

in) 

1.783 

cm 

(0.702 

in) 

1.783 

cm 

(0.702 

in) 

2.54 

cm 

(1.00 

in) 

0.76 

cm 

<0.30 

in) 

0.767 

cm 

(0.302 

in) 

0762 

cm 

(0.300 

in) 

0.777 

cm 

(0.306 

in) 

0.767 

cm 

(0.302 

in) 

0.76? 

cm 

(0.302 

in) 

0762 

cm 

(0.300 

in) 

1.27 

cm 

(0.50 

in) 

1.278 

cm 

(0.503 

in) 

1.267 

cm 

(0.499 

in) 

1.270 

cm 

(0.500 

in) 

1.267 

cm 

(0.499 

in) 

1.295 

cm 

(0.510 

in) 

1.265 

cm 

(0.498 

in) 

1.78 

cm 

(0.70 

in) 

1.786 

cm 

(0.703 

in) 

1778 

cm 

(0.700 

in) 

1.783 

cm 

(0  702 

in) 

1.775 

cm 

<0.699 

in) 

1.773 

cm 

(0.698 

in) 

1.770 

cm 

(0.697 

in) 

3.81 

cm 

(1.50 

in) 

0.76 

cm 

(0.30 

in) 

0.775 

cm 

(0.305 

in) 

0765 

cm 

(0.301 

in) 

0.775 

cm 

(0.305 

in) 

0.765 

cm 

(0301 

in) 

0.772 

cm 

(0.304 

in) 

0767 

cm 

(0.302 

in) 

1.27 

cm 

(0.50 

in) 

1.280 

cm 

(0.504 

in) 

1.275 

cm 

(0.502 

in) 

1.283 

cm 

(0.505 

in) 

1.273 

cm 

(0.501 

in) 

1.270 

cm 

(0.500 

in) 

1.260 

cm 

(0.496 

in) 

1.78 

cm 

(0.70 

in) 

1.786 

cm 

(0.703 

in) 

1786 

cm 

(0700 

in) 

1.778 

cm 

(0.700 

in) 

1773 

cm 

(0.698 

in) 

1.793 

cm 

(0.706 

in) 

1  781 

cm 

(0.701 

in) 

5.08 

cm 

(2.00 

in) 

0  76 

cm 

(0.30 

in) 

0.777 

cm 

(0.306 

in) 

0.762 

cm 

(0.300 

in) 

0.765 

cm 

(0.301 

in) 

0.757 

cm 

(0.298 

in) 

1  27 

cm 

(0.50 

in) 

1.273 

cm 

(0  501 

in) 

1.260 

cm 

<0.496 

in) 

1.275 

cm 

(0.502 

in) 

1.265 

cm 

<0.498 

in) 

1.283 

cm 

(0.505 

in) 

1.273 

cm 

(0.501 

in) 

1.78 

cm 

(0.70 

in) 

1.796 

cm 

(0707 

in) 

1765 

cm 

<0.695 

in) 

1.791 

cm 

(0.705 

in) 

1778 

cm 

(0700 

in) 

1  770 

cm 

<0.697 

in) 

1763 

cm 

(0  694 

in) 

Diameter  Imper 
faction.  J  ^ 
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